Modulation of Luttinger liquid exponents in multi-walled carbon nanotubes 
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We develop in this paper a theoretical framework that applies to the intermediate regime between 
the Coulomb blockade and the Luttinger liquid behavior in multi-walled carbon nanotubes. Our 
main goal is to confront the experimental observations of transport properties, under conditions 
in which the thermal energy is comparable to the spacing between the single-particle levels. For 
this purpose we have devised a many-body approach to the one-dimensional electron system, incor- 
porating the effects of a discrete spectrum. We show that, in the crossover regime, the tunneling 
conductance follows a power-law behavior as a function of the temperature, with an exponent that 
oscillates with the gate voltage as observed in the experiments. Also in agreement with the experi- 
mental observations, a distinctive feature of our approach is the existence of an inflection point in 
the log-log plots of the conductance vs temperature, at gate voltages corresponding to peaks in the 
oscillation of the exponent. Moreover, we evaluate the effects of a transverse magnetic field on the 
transport properties of the multi-walled nanotubes. For fields of the order of 4 T, we find changes 
in the band structure that may be already significant for the outer shells, leading to an appreciable 
variation in the power-law behavior of the conductance. We then foresee the appearance of sensible 
modulations in the exponent of the conductance for higher magnetic fields, as the different subbands 
are shifted towards the development of flat Landau levels. 

PACS numbers: 73.63.Fg, 73.22.-f, 73.23.-b 



I. INTRODUCTION 

Recent progresses in nanotechnology has revealed 
many novel transport phenomena in mesoscopic low- 
dimensional structures. Carbon nanotubes (CNs) were 
discovered by Iijima in 1991 1 -, as a by-product of carbon 
fullerene production. Thus, since then, a new field of 
research in mesoscopic physics has opened^. CNs are ba- 
sically rolled up sheets of graphite (hexagonal networks of 
carbon atoms) forming tubes that are only nanometers 
in diameter. The electronic properties of CNs depend 
on their diameter and chiral angle (helicity) parameter- 
ized by a roll-up (wrapping) vector (n, m)&. Hence it 
follows that some nanotubes are metallic with high elec- 
trical conductivity, while others are semiconducting with 
relatively low band gaps. They may also display dif- 
ferent electronic properties depending on whether they 
are single-walled carbon nanotubes (SWNTs) or multi- 
walled carbon nanotubes (MWNTs). MWNTs are typ- 
ically made of several (typically 10) concentrically ar- 
ranged graphene sheets with a radius of about 5 nm and 
lengths in the range of 1 — 100 /im. 

Several electrical transport experiments have been 
done by using CNs. In these experiments the effects 
of the Coulomb interaction are quite relevant, and the 
way they manifest depend mainly on the size of the sys- 
tem, the temperature, and the quality of the contacts 
used in experiments. Thus, MWNTs show for instance 
different transport regimes from ballistic transport (con- 
ductance quantization^, Luttinger liquid-like^) to dif- 
fusive transport and also strong localization&iS. 

When the contacts are not highly transparent, the 
measurements of the conductance and the differential 



conductivity reflect the strong Coulomb repulsion in 
CNsii*i£. At low temperature, T < 1 K, the zero-bias 
conductance is strongly influenced by quantum mechan- 
ical effects and shows regular oscillations corresponding 
to the so-called single electron tunneling (SET), because 
it arises from the electrons going one at a time through 
the devic o 13 : 14 . The voltage between two peaks is deter- 
mined by the energy for adding an electron to the device 
between the potential barriers. This addition energy de- 
pends on the level spacing, due to the momentum quanti- 
zation, and on the electron-electron interaction**. This is 
the Coulomb blockade (CB) regime where the tunneling 
transparency of a barrier vanishes owing to the electron- 
electron interaction. The differential conductivity of in- 
dividual MWNTs has been measured in Ref. Ha. where 
CB and energy level quantization have been observed. 

When the temperature increases, so that the thermal 
energy is much larger than the level spacing, the trans- 
port measurements reflect instead the many-body prop- 
erties of the system. In fact, the electron-electron interac- 
tion in a one-dimensional (ID) system is expected to lead 
to the formation of a Luttinger liquid (LL) with proper- 
ties very different from those of the non-interacting Fermi 
gasii. The LL state is characteristic of ID clean (ballis- 
tic) conductors, which have unusual electronic properties 
that cannot be explained by Fermi liquid theory. 

Landau quasiparticles are unstable in ID systems of 
interacting electrons, in which the low-energy excitations 
take the form of plasmons (collective electron-hole pair 
modes) that are the stable excitations: this is known as 
the breakdown of Fermi liquid in iDiSiiiLSl. The LL state 
has two main features: 

i) the power-law dependence of physical quantities, 
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such as the tunneling density of states (TDOS), as a func- 
tion of energy or temperature; 

ii) the spin-charge separation: an additional electron 
in the LL decays into decoupled spin and charge wave 
packets, with different velocities for charge and spin. 

Characteristic experimental signatures support the as- 
sumed LL behaviour of CNs 21 - 22 - 23 - 24 . In fact the power- 
law dependence of physical observables follows from the 
behavior of the TDOS as a function of the energy or 
the temperature: the tunneling conductance G reflects 
the power law dependence of the TDOS in a small-bias 
experiment^ 

G = dl/dV cx r (1) 

for eVb -C fcsT, Vb being bias voltage. Such kind of 
signatures of LL behavior have been observed in single- 
walled as well as in multi- walled nanotubes. In MWNTs 
the observed values of a lie between 0.04 and 0.3, and 
results in agreement with the experimental values have 
been obtained in previous theoretical work32S*SLSSi22. 

In this paper we introduce a theoretical description 
of the crossover as the temperature decreases from the 
LL regime to the CB regime. We do this in order to con- 
front the experimental data reported in a recent letter by 
Kanda et alm&, in which the intermediate regime has been 
explored measuring the zero-bias conductance at temper- 
atures where the thermal energy becomes comparable to 
the level spacing in the discrete single-particle spectrum. 
In Ref. |3(J the authors have reported a systematic study 
of the gate voltage dependence of the LL-like behaviour 
in MWNTs, showing the dependence of the exponent a 
on gate voltage, with values of a ranging from 0.05 to 
0.35. The main results of Ref. [3(] are as follows: 

i) the gate-voltage (V^) dependence of the exponent a 
below 30 K exhibits periodic oscillations; the character- 
istic V g scale for a variation, AV g , is around IV; 

ii) changes in the exponent a are observed in the plots 
of the conductance at an inflection temperature T* ~ 30 
K, for values of V g corresponding to peaks of a. 

iii) the exponent a depends significantly on the trans- 
verse magnetic field acting on the CN. a is reduced from 
a value of 0.34, at a peak in the a oscillations, to a value 
0.11, for a magnetic field B = 4 T. In most cases, the 
exponent a becomes smaller for higher magnetic fields. 

In their letter, Kanda et al. do not find plausible the 
explanation of the mentioned features starting from a LL 
description of the MWNTs. In this respect, Egger has 
considered in Ref. OH] a model for MWNTs composed of 
a number Nsh of ballistic metallic shells. The author 
has discussed there a low-energy theory for the MWNTs 
including long-ranged Coulomb interactions and internal 
screening effects. The theory may be also extended to in- 
clude the effect of a variable number of conducting modes 
modulated by the doping level. However, Kanda et al. 
rule out the possibility that the change in the number 
of subbands at the Fermi level may be at the origin of 
the features observed in their experiment, as long as the 



subband spacing is too large to be consistent with the 
period of the oscillations. 



Then, Kanda et al. turn to a different kind of the- 
ory that considers the MWNT as a diffusive conduc- 
tor. Egger and Gogolin2S have calculated the TDOS of 
doped MWNTs including disorder (with mean free path I 
smaller than the radius R) and electron-electron interac- 
tions. MWNTs may display an effective and nonconven- 
tional CB arising from tunneling into a strongly interact- 
ing disordered metal, leading to LL-like zero-bias anoma- 
lies: the exponent becomes a = (R/2itTiDvq) log(l + 
vqU), where D is the diffusion constant, v = N s /2ithvp 
is the noninteracting density of states depending on the 
number of subbands N s and the Fermi velocity Vp, and 
Uq is an effectively short-ranged ID interaction. Sub- 
stituting these parameters by pertinent values yields an 
estimate a ~ R/N s l, that is near the experimentally ob- 
served values. 



Kanda and coworkers conclude that this second the- 
ory better explains the experimental results, by assum- 
ing that the mean-free path / may fluctuate with the gate 
voltage. This is based on the theoretical work by Choi et 
al~, that have studied the effects of single defects on the 
local density of states via resonant backscattering. How- 
ever, the argument of Kanda et al. has not addressed 
the question of how a random distribution of defects may 
produce the oscillations observed in the experiment. We 
believe otherwise that the depedence of the a exponent 
on the gate voltage is in correspondence with a definite 
periodic structure of the single-particle density of states. 



In the next section, we put forward the idea that there 
exists a clear relation between the periodic oscillations in 
the a exponent and the quantization of the energy lev- 
els in the MWNTA Building on this idea, in section 
III we introduce an energy distribution function which 
takes into account the effects of the quantization at low 
temperature and can be compared with the Fermi-Dirac 
distribution in the limit of high temperature. Next, in 
section IV we develop a theoretical approach based on 
many-body theory, in order to deal with the low-energy 
effects of the the long-range Coulomb interaction in the 
ID electron system. We introduce in that scheme the dis- 
cussed energy distribution function, for the sake of evalu- 
ating the effects of the discrete spectrum on the exponent 
a. The results of this approach are reported in section 
V, where we carry out the comparison between our the- 
oretical predictions and the experiments. The effects of 
a transverse magnetic field are considered in section VI. 
In section VII we highlight how the distinctive features 
of our approach match naturally with the experimental 
observations, and we close with some concluding remarks 
and perspectives in section VIII. 
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II. ANALYSIS OF THE EXPERIMENTAL 
OSCILLATIONS 

Here we discuss the values of relevant parameters, in 
order to compare later the experimental results with the 
theoretical predictions. 

The transverse quantization energy of the nanotube, 
which corresponds to the distance between two nearest 
subbands, is given by AE± Tlvf/R, where the Fermi 
velocity is vf = 8 x 10 5 m/s and we take a nanotube 
radius R w 5 nm. On the other hand, the longitudinal 
quantization energy reads Ael = hvp/L, where in the 
experiments L w 4 -j- 6/im (here we are considering the 
entire length of the MWNT, L « 5.7(4.7) fj,m, rather 
than the length of its portion between the electrodes, 
d 0.7(1.0) /zm). The ratio turns out to be w 

^« 10 -100. 

In comparing the theoretical estimates with the ex- 
periment, there appears a conversion factor given by the 
coupling constant 7 between the gate electrode and the 
MWNT. This factor is approximately 7 = C/Cq, where 
C is the capacitance between the metal electrodes and 
the MWNT while Cq is the capacitance between the gate 
electrode and the MWNT. In the experiments by Kanda 
et al£Mk C « IfF and C w laF so that we obtain 
7 « 1 x 10 3 . 

Using these parameters Kanda and coworkers confute 
the above mentioned model based on the LL theory. In 
fact, the transverse quantization energy is about 0.1 eV, 
so that the typical value AV g = -fAE± for the change in 
the number of subbands N s at the Fermi level is about 
10 2 V. This value is much larger than the period for the 
gate voltage in the observed oscillations, which corre- 
sponds to some Volts. This theory would also predict 
a stepwise change in a as a function of V g , which is not 
in agreement with the oscillatory behavior measured. 

In the experiment, the measured value of the period 
AV g is around 3 V at a temperature T below 30 K. 
As shown by the cited authors in a previous paper—, 
where they reported the observation of the CB effect in 
MWNT, a similar value of AV g corresponds to a period 
(a sequence of 4 peaks or diamonds) in the zero bias con- 
ductance at very low temperature T ~ 30 mK. 

Following the usual formalismi* we can calculate the 
period as AV^ = ^(vph/L + 4U) 16 . If we neglect the 
electron-electron interaction U, we can assume AV g ~ 
"fVph/L, in order to compare the number of oscillations 
in a definite range of the gate voltage (— 10V to 10V 
in the experiment) observed for two different nanotubes 
(i.e. samples A and B in Fig. 1, respectively). Thus, 
we expect that the number of oscillations depends on the 
length as 
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From the experimental data (see Fig. 1) we observe re- 
spectively 10 and 8 oscillations, corresponding to an en- 



FIG. 1: A figure summarizing the gate-vo ltag e dependence of 
a in sample A (La = 5.7/xm top) of Ref. l3fl calculated from 
the data below 30° -K\ A similar a — V g plot is obtained for 
sample B (La = 4.7/im bottom). 



ergy spacing Ae w 2—2.5 meV, in perfect agreement with 
our estimate. This agreement suggests a theory based on 
a microscopic description of the single-particle spectrum, 
in order to explain the oscillations. Here the quantized 
energy levels play a central role. This theory has to be 
valid in the intermediate regime between very low tem- 
perature, where the peaks of the CB are measured by 
increasing V g , and room temperature, where the LL be- 
haviour prevails. To this aim, we will rely on the ID 
many-body approach, incorporating at the same time the 
discrete character of the spectrum which arises from the 
finite geometry of the system. 



III. SINGLE-PARTICLE PROPERTIES 

In the introduction of microscopic features in the 
many-body formalism, a central role is played by the en- 
ergy distribution function (EDF) rt(e;£p). In the many- 
body approach to electron liquids, the function n(e,SF) 
corresponds to the Fermi-Dirac distribution, which is 
usually replaced by the Heaveside step function — ef) 
in order to simplify the calculations. 

At very low temperatures, as the ones that are needed 
for the CB regime, the EDF has to reflect the discrete- 
ness of the energy levels. Thus, some pattern of peaks, 
corresponding to the quantized energy levels, has to be 
present. Hence, at low temperatures, where the micro- 
scopic level quantization is manifest, peaks arise with a 
shape smoothed by the growth of the temperature. The 
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spacing Ae of the peaks depends on the nanotube length 
L 



Ae 



V F h 



This value yields a temperature T c « As /kg, above 
which we expect that the oscillations will be smoothed. 

The thermal energy fcgT c governs then the crossover 
between two different regimes. From our analysis, it is 
clear that two scales for the behavior as a function of 
temperature characterize the system, corresponding to 
the longitudinal and the transverse quantization energies, 
although the second one turns out to be about 10 -j- 100 
times larger than ksT c . 

We can model the n-th peak of the EDF by a normal- 
ized function /(e,n), with a dependence on temperature 
consistent with the line shape for a thermally broadened 
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4fcTcosh^ 



' e-nAe ' 
v 2kT , 



(2) 



where nAe stands for the position of the n-th energy 
level. The energy distribution function is cutoff by the 
Fermi energy s F , as it is given by 



i{e; e F ) = ^ n ) °( e F ~ £ )- 



(3) 



n=0 



Alternatively we could introduce a smooth cutoff at the 
Fermi energy, taking into account thermal effects 



n(e;e F ) 



E 

n=0 



1 



1 



,(e-s F )/kT ' 



(4) 



Anyway the sharp cutoff given by 9{e F — e) is more conve- 
nient, in order to proceed with analytical computations. 
On the other hand, we have checked that the presence of 
the smooth cutoff does not change appreciably the whole 
scenario. 

At high temperatures, we can neglect the level spacing 
and calculate the EDF by integrating 



n(e\e F ) 
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dyf(e,y) 
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1 -|- e (e-e F +Ae/2)/kT ' 



(5) 



The latter formula gives the Fermi-Dirac distribution, 
which is usually approximated with the step function. 
However, the energy distribution function has in general 
a more structured shape, that depends on the ratio be- 
tween the level spacing Ae and the thermal energy kT as 
illustrated in Fig. 1. 

The structure of peaks in the energy distribution func- 
tion affects several objects in the many-body theory, like 
the one-loop polarization IT ) (q, cu q ). This object counts 
the number of particle-hole excitations that can be built 




FIG. 2: Plot of the energy distribution function n(e;ep) for 
different values of the temperature corresponding to kT / Ae = 
0.25, 0.6, 1.0 . The solid lines have been obtained from Eq. @ 
and the dashed lines by using the Fermi-Dirac distribution JSJ. 



with momentum q across the Fermi level. Thus, in the 
case of a model with linear dispersion e(k) ~ v F k, it 
becomes proportional at zero temperature to 

dk 6(q + k~ k F ) - 9{k - k F ) 
.2^ 



v F q 
1 f q dk 



1 



v F q J 2tt 2ttv f 



(6) 



In the approach based on the distribution function given 
by (J3J), the peaks give rise to a factor in the integration 
that reflects the discrete structure. Therefore Eq.© be- 
comes 

°° dk n(v F (q + k)\ e F ) — n(v F k; £ F ) 
-oo 27r v Fq 
1 r q dk 1 

= — / -AeJ2f(vFk,n)^- N T (e F ), (7) 

v F qJ 2tt ^ 2ttv f 

where at first order the function N F (e F ) can be taken in- 
dependent on q. Given the sharp cutoff, the N F function 
gets a simple form: 



/ dkAe f(v F k, n) 

J ° n=Q 



Ae ^ £f + VF q - nAe e F -nAe 
2^ Ej tanh ( 2kT } - tanh (^I^)] 



n=0 



2kT 



Ae ~ M , t/ e F -nAe,, 
W E^ 1 ~ tanh ( ^ )] = lN T {e F ) . (8) 
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Therefore it holds 
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FIG. 3: Plot of N T (e F ) for three different values of kT j Ae: 
kT/Ae =0.1 (dashed), =0.15 (solid), =0.3 (dotted-dashed) . 
e is in unit of Ae 



transferrin. It is known that the Coulomb interaction 
remains long-ranged in a ID electron systemSSiiS, and 
this property is also shared by the graphene sheet, which 
has a vanishing 2D density of states at the Fermi level. 
To obtain a sensible expression of the ID interaction, we 
may start from the representation of the l/|r| Coulomb 
potential in three spatial dimensions as the Fourier trans- 
form of the propagator 1 /k 2 

_l_ r d 3 k k . r j_ 

|r| J (27r)3 e k 2 [W) 

If the interaction is projected onto one spatial dimension, 
by integrating for instance the modes in the transverse 
dimensions, then the Fourier transform has the usual log- 
arithmic dependence on the longitudinal momentum^ 

i <*» J-iogCfc) (ii) 

\x\ J 2n 2n 



This function is represented in Fig. 3. As expected, 
N T (e) is periodic with period Ae. 

We introduced the function Nt depending on the 
Fermi energy ep. In this respect, the position of ep, 
tunable by the gate voltage V g , fixes the portion of the 
last peak that we have to integrate. It is clear that, from 
an experimental point of view, ep is in correspondence 
with the value of the gate voltage V g , so that it can be 
shifted in a continuous way, allowing for the exploration 
of the region between two peaks in Fig. 2. 

IV. MANY-BODY APPROACH 

In the CNs, the relevant interaction is given by 
the large Coulomb potential at small momentum- 



In our case, the scale of the longitudinal momentum k 
in the logarithm is dictated by the existence of a short- 
distance cutoff k c , from the small transverse size of the 
electron system. Thus, we end up with a representation 
of the Fourier transform V(k) of the ID Coulomb poten- 
tial 

V{k)*±-\og{k c /k) (12) 

where k c ~ l/R. 

Our interest focuses in describing the interaction 
among a number N s of subbands crossing the Fermi level 
in the conducting shell of a multi- walled nanotube. Thus, 
we write the hamiltonian for the electron fields $± off en- 
coding the modes of the different branches with linear 
dispersion about the Fermi points PF a , —PFa- 



H = dp Yj^taM {±V F p)i}±, 

"'-fee a =l 



In the above expression, p(p) stands for the charge den- 
sity made from the sum of all the electron densities for 
the fields ip± aa (p)- 

In the hamiltonian l|13|) . we have neglected interac- 
tions that lead to a change of chirality (±) of the fields, 
as they imply a large momentum-transfer of the order of 
2kp and are therefore subdominant with respect to the 
large Coulomb interaction at small momentum-transfer. 
Those so-called backscattering interactions are known to 
remain negligible down to extremely low energies. On 
the other hand, the model build from interactions be- 
tween currents with well-defined chirality has some inter- 



-{p)+e 2 [ ' dpp(p)V(p) p(-p) (13) 
J-fc c 

I 

esting properties. One of them is that, when performing 
the perturbative calculation of the polarization function 
n((7, ujq), the self-energy corrections cancel order by order 
against the vertex corrections. Thus, the exact result is 
just given by the bare polarization funcion. In our case, 
upon the introduction of the different subbands and the 
mentioned effects in the single-particle density of states, 
we have 

U(q,u q ) = -N s N T (e F ) " Fq \ (14) 
7T - v F q 2 

Moreover, there is further simplification of the dia- 
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FIG. 4: Self-consistent diagrammatic equation for the dressed 
Coulomb interaction. 



grammatics from the fact that all the fermion loops with 
more than two interaction lines vanish identically in ID. 
Thus, the dressed Coulomb interaction has to satisfy the 
self-consistent equation shown in Fig. 0] that amounts 
to the RPA sum of particle-hole contributions. 

We are mainly interested in obtaining the quasiparticle 
properties in the model governed by , which requires 
the computation of the electron self-energy H(q,u> q ). As 
long as the RPA provides the exact expression of the 
dressed Coulomb interaction, a sensible representation 
of the self-energy is given by the partial sum of the per- 
turbative expansion shown in Fig. |3J Correspondingly, 
we may express the electron self-energy in the form 

£(fc >Wfc ) =i[ kC ^ ^ %*G<°\k-q,w k - Uq ) 

J-k c 27r J-oo 27r 



l-V(q)R{q,u q y 



(15) 



where G^{k,u>k) stands for the free electron propaga- 
tor. For the evaluation of the scaling properties of the 
quasiparticle parameters, it is enough to obtain the loga- 
rithmic dependence of the function £((?, Uq) on the high- 
energy cutoff E c ~ v F k c . Thus, the representation in 
Fig. misses corrections to the fermion line G^ ' as well 
as corrections to the three-point vertex at the end of the 
line, but these do not introduce new divergences in the 
cutoff E c and can be safely disregarded in the determi- 
nation of the electron field scale. 



L(k,o>) = E 





n=l 



FIG. 5: Set of diagrams for the computation of the electron 
self-energy, where each particle-hole bubble represents the full 
polarization operator. 

We therefore concentrate on the computation of the 
contributions that diverge with the cutoff E c in l)15|l. 
Order by order in perturbation theory, we find that the 
same kind of logarithmic dependence appears in the lin- 
ear terms proportional to k and uj k . We can recast the 
whole perturbative series into the expression 



1 V{q )e 2 



2tt 3 / 2 
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71=0 
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x(cj k ± v F k) \og{E c 
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AN S N T y 
x{uj k ± v F k) log(.E c ) 



- 2 



(16) 



(17) 



where the effective coupling is given by g = 
2N s NtV (qo)e 2 /ttVf, in terms of a suitable average mo- 
mentum qy^. 

The frequency and momentum-dependence in Eq. I|17|) 
shows that such self-energy corrections represent a renor- 
malization of the electron field operator. Then it is con- 
venient to define a renormalized electron field \P r related 
in each case to the bare electron field >F through a scale 
factor Z 1 / 2 (E C ): 



= Z 1 ' 2 {E C )^! 



(18) 



The renormalization factor Z(E C ) must be fixed so that 
it makes finite the renormalized electron Green function 
G(k,u>k)' This is given by 
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G 
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- E(jfe,wfc) 



G(°) 

Z-\E c )(Lu k ±v F k) 
+Z-\E c ){u k ± v F k)a{g) \og(E c ) 



where 



1 



AN S N T 



(19) 



(20) 



The requirement of cutoff- independence of the electron 
Green function leads to absorb the logarithmic depen- 
dence on E c into the renormalization factor Z(E C ). Thus 
we get the scaling equation 



E c —logZ(E c ) = a(g). 



(21) 



We observe that Z(E C ) follows a power-law behavior as 
a function of the energy cutoff E c . We may identify the 
factor Z(E C ) with the weight of the electron quasipar- 
ticles, which vanishes as the Fermi level is approached 
in the limit E c — > 0. We capture in this way one of 
the most genuine features of the many-body theory in 
ID, namely the absence of low-energy excitations with 
fermionic character, which defines the Luttinger liquid 
class of ID electron liquids. 



V. TUNNELING DENSITY OF STATES 

The power-law behavior of the quasiparticle weight 
Z(E C ) translates into a similar dependence on energy of 
the tunneling density of states n(s). By trading the de- 
pendence on the high-energy cutoff by the energy mea- 
sured with respect to the Fermi level, we obtain the scal- 
ing behavior 



n(e) - N s Z{s) - e a{g) 



(22) 
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Usually the behavior of the TDOS is obtained by measur- 
ing its dependence on temperature or bias voltage. The 
exponent in (|22|) is then to be compared with the mea- 
surements of power-law behavior observed in the plots 
of conductance versus temperature, or in the differential 
conductivity. 

Regarding the experiments in Ref. we are in dispo- 
sition to explain the oscillations observed in the exponent 
of the conductance, in a temperature regime correspond- 
ing to a thermal energy of the order of the level spacing 
Ae. The point is that, as we shift the Fermi level along 
the structure of single-particle peaks in the spectrum, the 
a exponent has to reflect the modulation of the periodic 
factor Nt(s). As mentioned before, the period of the os- 
cillations is in agreement with the change in gate voltage 
required for a shift Ae of the Fermi level. We comment 
now on several features regarding the amplitude of the 
oscillations. 

It becomes very convenient to plot the expression of 
the exponent a(g) as a function of the variable N s Nt- 
This depends also on the effective coupling V(qo)e 2 /vf of 
the Coulomb interaction. We have estimated this quan- 
tity by using a value of qo consistent with the nanotube 
length, and taking into account the screening of the in- 
teraction, mainly by other shells in the nanotube. Thus, 
sensible values of V(qo)e 2 /vf turn out to be of the order 
of ~ 10. The exponent a has been represented in Fig. 
as a function of N s Nt, for two different choices of the 
effective coupling. We may assume the typical situation 
in which the MWNTs are hole-doped, with a number of 
subbands crossing the Fermi level iV s « 6 - 10 41 . We 
observe that, for such an estimate of iV s , the modulation 
of the Nt factor (with the amplitude reported in section 
III) accounts for oscillations of a reaching values as high 
as « 0.3 (for a large value of the effective coupling) and as 
low as ~ 0.05 (for a small value of g and a large number 
of subbands). 

It is actually quite plausible that the drift observed in 
general in the average value of the a exponent (towards 
larger values for higher gate voltages) may be the con- 
sequence of a change in the Fermi velocity as the Fermi 
level approaches the top of the lowest partially filled sub- 
band. We note that the height reached by the oscillations 
in the a exponent is different in the two plots of experi- 
mental measures shown in Figs. 2(c) and 2(d) in Ref. I30L 
This may be attributed to a different relative position of 
the Fermi level between two consecutive subbands in the 
nanotube bandstructure. We have modeled the proxim- 
ity to the point with divergent density of states at the 
top St of a given subband by allowing for a dependence 
of the Fermi velocity vf{zf) ~ V £ t — £f- When this is 
taken into account in conjunction with the modulation of 
the periodic function Nt(ef), we obtain a picture for the 
dependence of the a exponent on gate voltage which is in 
good qualitative agreement with the oscillations reported 
in Ref. |3(], as shown in Fig. [7| 

We observe the similarity of the plot in Fig. [7| with 
the experimental curves in Fig. 1 for the exponent over 
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FIG. 6: Plot of the exponent a as a function of the product of 
the number of subbands N s and the weight Nt in the single- 
particle spectrum. 




FIG. 7: Plot of the exponent a as a function of the Fermi 
energy ef (in units of Ae) for different values of temperature 
kT/Ae = 0.1 (dashed), 0.15 (solid), 0.3 (dotted-dashed) . 



the 20 V range in the gate voltage. The oscillations are 
significant for low values of the thermal energy fcsT, as 
compared to Ae. For higher temperatures the thermal 
fluctuations erase the structure of superposed peaks in 
(|3*|) , what leads in turn to the disappearance of the oscil- 
lations in the a exponent. 

An important point reinforcing the above interpreta- 
tion of the experimental measures is that the depeden- 
dence on the variable N s Nt(£f) shown in Fig. trans- 
lates naturally into an asymmetric behavior of the ex- 
ponent a, when the oscillations approach the peak of 
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the curve near a ss 0.3. If we consider for instance a 
doped nanotube with N s = 10 and an effective coupling 
V(qo)e 2 jvp w 12, we observe that the exponent may run 
from a value a « 0.18 for Nt = 1 to a value a » 0.3 
for reasonably small Nt- However, it is clear that the 
other side of the oscillation does not have the same mag- 
nitude. This is consistent with the fact that inflection 
points in the plots of the conductance versus tempera- 
ture have been observed for peak values in the exponent 
oscillations. In our model, the values corresponding to 
a dip for Nt > 1 do not deviate much from the high- 
temperature behavior recovered for Nt ~ 1. It is also 
remarkable that the inflection points have been found 
experimentally to correspond to the highest oscillations 
for large gate voltage. Again, this is consistent with our 
theoretical interpretation, in which the lower values of 
the exponent a correspond to smaller values of the effec- 
tive coupling that do not support large variations upon 
changes in Nt, as shown from the lower curve in Fig. [5] 

VI. MAGNETIC FIELD EFFECTS 








-0.02 -0.01 0.01 0.02 



momentum 

FIG. 8: Energy dispersion relation of a (280,0) zig-zag nan- 
otube with B — (solid lines) and B — 4T (dashed lines). 
The position of the Fermi level for a hole-doped nanotube is 
also indicated. Energy is in eV (we used a hopping param- 
eter t = 2.5eV between nearest neighbor carbon atoms) and 
momentum is in A" 1 ). 



The experiments reported in Ref. |53 have also shown 
that the exponent in the power-law behavior of the zero- 
bias conductance may depend on the applie d tran sverse 
magnetic field. The magnetic length I — ^Jc/eB corre- 
sponding to the field strength used in the experiments 
(B = 4 T) is anyhow larger than the radius of the outer 
shells in the MWNTs. Then, the resulting band struc- 
ture must be far from the development of Landau levels, 
although it must already lead to the features producing 
the variation observed in the a exponent. 

To study the influence of the magnetic field, we have 
analyzed the change that it may produce in the band 
structure of nanotubes with different chiralities. Previ- 
ous investigations have addressed the formation of Lan- 
dau subbands in carbon nanotubes upon application of 
a transverse magnetic field. Here we have undertaken 
a similar computational approach, but applied to nan- 
otubes of large radius R ~ 10 nm, corresponding to ac- 
tive shells in multi- walled nanotubes. We have thus re- 
sorted to a tight-binding calculation of the hamiltonian 
in the carbon lattice, with the usual prescription of cor- 
recting the transfer integral by appropriate phase factors 

expfi^jf Adl^J (23) 

depending on the vector potential A between nearest- 
neighbor sites r and r'. The vector potential has been 
chosen in particular with the appropriate functional de- 
pendence at the nanotube surface, as described in Ref. 

m 

The results of the numerical diagonalization of the 
tight-binding hamiltonian for a magnetic field B = 4 T 
are shown in Fig. [H] The band structure corresponds to 
the particular case of a zig-zag nanotube, but it is quite 



representative of the changes produced in chiral nan- 
otubes upon application of the magnetic field. There is 
an incipient tendency of the subbands to flatten, an effect 
which is precursor of the formation of Landau subbands 
at much higher magnetic field. The two low-energy sub- 
bands spanning the gap for semiconducting nanotubes 
tend to approach, as the lowest Landau level is always 
at zero energy in the carbon nanotubes. The next sub- 
bands move in the opposite direction, towards a very slow 
reaccomodation along the direction of the higher Landau 
levels. 

In this perspective, there is a significant modification 
in the low-energy spectrum whenever the number of sub- 
bands crossing the Fermi level changes, as a consequence 
of switching on the magnetic field. This happens in par- 
ticular when the Fermi level for B = is slightly below 
the top of one of the subbands which is pushed down- 
wards under the action of the magnetic field. Recalling 
the discussion of the preceding section, this is the in- 
stance that corresponds to the higher gate voltages in 
Fig. [7| where the value of the exponent a w 0.3 is ap- 
proached. In these circumstances, a couple of degenerate 
subbands with large density of states are lost at the Fermi 
level upon application of the magnetic field. Although 
the number N s of subbands contributing to the conduc- 
tion decreases, the main effect influencing the computa- 
tion of the exponent a comes from the larger value of vf 
in the remaining (outer) degenerate subbands. Taking a 
value of the Fermi velocity in correspondence with the 
Fermi level shown in Fig. [3] we find that the exponent is 
reduced down to a value a ~ 0.1, in agreement with the 
experimental observation reported in Ref. l30l 

Although no measures have been reported over a given 
range of gate voltage in the presence of magnetic field, it 
is very instructive to apply Eq. (|20|) to find the results 
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FIG. 9: Plot of the exponent a as a function of the Fermi 
energy ef (in units of Ae) for B = (solid line) and B — 4T 
(dashed line). Here the temperature is kT/Ae = 0.15. 

from the change in the band structure shown in Fig. [5] 
The gate- voltage dependence obtained in this way for the 
a exponent at B = 4 T is represented in Fig. |5J We ob- 
serve that the plot of the exponent is now below the curve 
corresponding to the oscillations at B = 0, which is in 
agreement with the general statement made in Ref. l30l 
However, the comparison with the experimental obser- 
vation reported at a dip in the oscillations (a ~ 0.06 at 
zero magnetic field) shows that our prediction falls short 
to account for the measured reduction. The experimental 
result reported at B = 4 T, a w 0.005, implies anyhow 
a too large effect on transport properties, which is diffi- 
cult to reconcile with the change in the band structure 
for such a relatively small magnetic field. In this respect, 
we point out that the evidence for a well-defined contin- 
uous plot is not so clear for that particular measure of 
the conductance as in other cases^. More experimental 
input will be needed, over a complete range of gate volt- 
ages, in order to confront the theoretical predictions with 
a more precise picture of the magnetic field effects on the 
transport properties. 



VII. DISCUSSION 

Now we discuss the main distinctive features of our 
theoretical approach, which lead to very good agreement 
with the experimental observations. 

I) The novelty of our framework reflects itself in the 
prediction of the periodic dependence of the exponent a 
on the Fermi energy of the system, due to the oscillatory 
behavior of the factor Nt(ef) as the Fermi level sweeps 
the structure of single-particle levels. We estimated the 
value of Ae in section II and inferred the variation needed 
in the gate voltage V g to produce that shift in energy 
from the capacitive coupling between the gate and the 
MWNT. We thus estimated a value of the period for V g 



of the order of ~ 2 V, in agreement with the experi- 
mental results. As we have shown in Fig. 7, our model 
can explain the oscillations in the a exponent, leading to 
plots that are in good qualitative agreement with the ex- 
perimental results showing the dependence on the gate 
voltage. The precise form of these experimental obser- 
vations may change from one sample to another, but it 
remains clear that the oscillations with periodicity AV g 
are in correspondence with the level spacing due to the 
quantization of single-particle levels. 

II) The oscillations in the experimental plot of a vs. V g 
were observed for values of ksT lower than a threshold 
which is consistent with our estimate for the level spacing 
Ae. Furthermore, our approach explains naturally the 
existence of an inflection point T* in the log-log plot of 
the conductance versus the temperature, for values of V g 
corresponding to peaks in the plot of a. As reported in 
Ref. l30t the conductance keeps following an approximate 
power- law behavior below the temperature T* , but with 
an exponent significantly enhanced with respect to the 
high-temperature value. In our model, this is a natural 
consequence of the asymmetry in the range of variation 
of a in the upper plot of Fig. [|J1 when the mean value of 
the product N s Nt is in the range 6 — 10. Peak values of 
a as large as « 0.3 can be reached at low temperatures 
(which means Nt < 1). However, minimum values of 
a are reached at the other side of the oscillation, for 
Nt > 1. In that case, low temperature does not mean 
a significant change of the exponent with respect to the 
mean, high-temperature value. 

Thus, it can be shown that our approach provides a 
sensible description of the power-law dependence and the 
inflection point in the plot of the conductance as a func- 
tion of temperature. In our framework, we may introduce 
the temperature as the relevant variable by trading the 
scale depedence on the cutoff E c by scale dependence on 
the thermal energy fcgT. Then we can apply renormal- 
ization group equations like H21[) . but now taking into 
account that the right-hand-side a(g) depends on tem- 
perature through the factor Nt- We have represented in 
Fig. ^3 the behavior of the conductance G that results 
from the temperature dependence of the renormalization 
factor Z(fcsT). We observe that an inflection point in 
the log-log plots appears only for peak values of the a 
exponent, with values of the crossover temperature T* 
corresponding to ksT* ~ As^. Our description is also 
consistent with the experimental fact that such inflection 
points arise for gate voltages giving the largest values of 
a. In our approach, the smaller values of a correspond 
in general to oscillations along the lower curve in Fig. 
[51 which naturally lends less sensitivity to the exponent 
upon changes in the temperature variable. 



VIII. CONCLUDING REMARKS 

In this paper we have developed a framework that 
applies to an intermediate regime between that of the 
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FIG. 10: Log-log plots of the conductance as a function of 
the temperature. In the left plot the Fermi energy is fixed at 
sf ~ 3.7Ae, i.e. at the last peak shown in Fig. Q I n the 
right plot the Fermi energy is fixed at ef ~ — 3.7Ae, i.e. at 
the first minimum shown in Fig. |7| 

Coulomb blockade and the Luttinger liquid behavior. 
Our approach accounts for the main experimental fea- 
tures in the zero-bias conductance reported in Ref. I30L 
This represents a significant progress with respect to the 
arguments presented in that reference to interpret the ex- 
perimental observations. There it was actually claimed 
that a nonconventional Coulomb blockade, arising from 
tunneling into a strongly interacting disordered metal, 
could be at the origin of the exponent oscillations. How- 
ever, it was also noted that such a theory does not predict 
power-law behavior above the inflection temperature T* . 
The authors of Ref. |3fJ concluded that "detailed theoret- 
ical and experimental investigations would be required 
for quantitative consideration" . 

We believe actually that, regarding transport measure- 
ments, the role of disorder is different than that at- 
tributed in Ref. |33 to interpret the experimental re- 
sults. The nonconventional Coulomb blockade predicts 
a power-law behavior of the conductante, with an expo- 
nent given in terms of the mean free path ibyaw R/N s l. 
Within this theory, the oscillatory behavior of the expo- 
nent was related to the possible modulation of the mean 
free path coming from disorder. This was supported by 
recalling the theoretical work of Choi et alm& about the 
influence of different forms of disorder on the local density 
of states. A closer look at this study shows, however, that 
it is highly unlikely that the combined effect of different 
defects may produce the desired oscillatory pattern. 

From the results of Choi et al. in Ref. 1331 we note 
that most part of the impurities and defects of the carbon 
lattice have an effect on the density of states extending 



over an energy range much larger than 1 meV. This is 
the energy scale consistent with the observed oscillations 
in the a exponent. Thus, from that source of disorder 
we can only expect a more or less uniform effect on the 
range of gate voltage shown in the plots of Fig. 1. On 
the other hand, the effect of vacancies in the lattice is 
to produce much narrower peaks in the density of states, 
which now may give rise to features discernible at the 
meV scale. We believe that this may be the main origin 
of the irregularities observed in the oscillations of the a 
exponent. 

Finally, we have also shown that our framework is con- 
sistent with the experimental results reported in Ref. 
about the magnetic field effects in the transport prop- 
erties. Using a tight-binding approximation adapted to 
include the action of a transverse magnetic field, we have 
computed the change in the band structure of the thick 
shells in multi-walled nanotubes. In this case, the radius 
of the shells becomes comparable to the magnetic length 
already for a magnetic field of 4 T. Then, we have seen 
that the low-energy subbands show an incipient tendency 
to flatten along Landau levels. This may induce in gen- 
eral further modulation of the exponent a, depending on 
the strength of the magnetic field. 

We have seen that, for high values of a at B = (about 
~ 0.3), our approach predicts a sensible reduction of the 
exponent, which is in agreement with the experimental 
observation. As long as the density of states at the Fermi 
level is reduced upon application of the magnetic field, we 
find a picture that is consistent with the statement that, 
in most cases, a is smaller for higher magnetic fields^. 
However, a clarification from the experimental point of 
view is required, in order to know the exceptions to that 
rule. In our framework, we may envisage some instances 
in which the application of the magnetic field may in- 
duce the proximity of the Fermi level to the top of a 
given subband, with the consequent increase in the den- 
sity of states (and the effective e-e interaction). We find 
in this way that sufficiently large variations in the mag- 
netic field may produce a modulation in the exponent a, 
rather than a monotonous trend. Our model provides 
anyhow definite predictions that are susceptible of being 
confronted experimentally. In this respect, more experi- 
mental input would be desirable to discern the effects of 
a transverse magnetic field, that may be quite significant 
for the outer shells in multi-walled nanotubes. 
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